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1 INTRODUCTION 



In the early Universe, just before recombination, fluctua- 
tions in the coupled baryon-photon fluid were subject to 
two competing effects: attractive gravity and repulsive pres- 
sure. These two effects are expected to produce a series 
of acoustic peaks - dubbed Baryon Acoustic Oscillations 

i BAOs) - in both the Cosmic Microwave Backgrou nd (CMB) 
Sunvaev fc ZeldovichI [l970l: [Peebles fc Yul ll970D and the 
matter power spectra (|Eisenstein fc Hul 19991 ). 

These features have been observed in the CMB 
temperature-temperature power spectrum using data from 
several years of the Wilkinson Microwave Aniso tropy Probe 
(WM AP) data (for the latest measurements see lNolta et al.l 
|2008| ). As galaxy surveys cover increasingly larger volumes, 
they too can probe the scales on which the BAOs are pre- 
dicted. Recently these h ave been observed in the Sloan Dig- 
ital Sky Survey (SDSS) (|Eisenstein et al.ll2005l : lHutsill2006al : 



ABSTRACT 

The Baryon Acoustic Oscillations (BAOs) or baryon wiggles which are present in the 
galaxy power spectrum at scales 100 — 150ft.^^Mpc are powerful features with which 
to constrain cosmology. The potential of these probes is such that these are now 
included as primary science goals in the planning of several future galaxy surveys. 
However, there is not a uniquely defined BAO Method in the literature but a range of 
implementations. We study the assumptions and cosmological performances of three 
different BAO methods: the full Fourier space power spectrum [P(/c)], the 'wiggles 
only' in Fourier space and the spherical harmonics power spectrum [C {£)] . We contrast 
the power of each method to constrain cosmology for two fiducial surveys taken from 
the Dark Energy Task Force (DETF) report and equivalent to future ground and space 
based spectroscopic surveys. We find that, depending on the assumptions used, the 
dark energy Figure of Merit (FoM) can change by up to a factor of 35 for a given fiducial 
model and survey. We compare our results with the DETF implementation and, discuss 
the robustness of each probe, by quantifying the dependence of the FoM with the 
wavenumber range. The more information used by a method, the higher its statistical 
performance, but the higher its sensitivity to systematics and implementations details. 

iPercival et al.ll07 a|) and the 2 degree Field Galaxy Survey 
(2dFGS) (|Percival et al.ll07 bh . 



The scale of the oscillations provides a standard 
ruler and has been used t o con s train dark ener g y pa- 
rameters (lEisenstein et al.l [2OO5I: lAmendola et al.| 2005 : 
Wang fc Mukherieell2006l:IWandl2006l : |Percival et al.llo7 al lb: 
Ichikawa fc Takahashil |2007^ as well as neutrinos masses 
CGopbar et al. 2006 !) and e ven alternative models of grav- 
itvlAlam fc Sahnill 2006: Pir es et al lliooel : ICuo et al.ll2006l : 
lYamamoto et al.l I2OO& : , WanJBoOSi ) . though current data 
does not cover enough cosmological volume to be constrain- 
ing without the help of external data sets. 

In the future, BAOs will be a fundamen tal tool for preci- 
sion cosmolog y (Amcn dola et al.| [2005: Pe acock et al.ll2006l : 
I Albrecht etTahl l2006i ) and there are many planned surveys 
which use BAOs as one of their primary science drivers. 
These include ground based surveys such as the Dark En- 
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ergy Survey llAnnis et al.l 2005|, DE Sl. the Large Synoptic 
Survey Telescope l|Zhan et al.BoOTl. LSST). the Wide -Field 
Multi-Object Spectrograph ('Glazcbrook et al.l l2005l . WF- 
MOS), WiggleZ (Gla zcbrook et al. 2007), the Bary on Os- 



cillations Spectroscopic Survey i Schlegel et alJlioOTl . BOSS 
or S PSS III), the Physics of the Accelerating Universe Sur- 
vey (|Benitej|2008l . P AU), the Hobby- Eberly Telescope Dark 
Energy Experiment (|Hill et al.ll2004l . HETDEX), as well as 
space-based surveys such as the Advanced Dark Energy 
Physics Telesc ope (ADEPT) the S pectroscopic All-sky Cos- 
mic Explorer (|Cimatti et al."2008', S PACE), the Dark UNi- 
verse Explorer (JRcfrcgicr & the DU NE collaboration! [iooj . 
DUNE) and the EUCLID projeclQ currently under study 
by the European Space Agency. Such surveys could probe 
the baryon oscillations in the gal axy di s tributi on as well 
as the distribution of clusters (|Hiitsil 2006bl). quasar s 
ISchlegel et al.ll2007l ) or even supernovae l|Zhan et al.ll2008l ). 

One of the attractive properties of using the BAOs as 
a cosmological probe is that they are considered to repre- 
sent a robust probe for extracting cosmological information. 
However there is not a unique BAO Method but a range of 
methods which differ in a number of ways. 

The first variable in the method is the statistics 
used in t he analysis: the meas urement can be done in 
real space llEisenstein et al.ll2005l . £,{r)), c onfiguration space 
(Lovordc ot al."2008', w{e)), Fourier space (|Seo fc EisensteinI 
12003 ; Amcndola ct a l, 2005 , P(k)) or in spherical harmonic 
space (jPolnev et al.l |2006| . C{1)). Even when using the 
same statistic, there are approaches which use different 
levels of information and are therefore subject to differ- 
ent systematics. For studies in Fourier space, some mea- 
sure t he full power spectrum information, including the 
BAOs (|Seo fc Eisensteinll2003l ). whereas others subtract the 
smoot h part of the spectrum and focu s only on the oscilla- 
tions l|Blake et al.ll2006l : ISeo fc Eisens tcin 2007) - the latter 
uses less information but may be more robust with respect 
to systematic uncertainties. 

In this paper we use the Fisher matrix approximation of 
the likelihood to contrast the information available in three 
of the methods described above: the Fourier space power 
spectrum [P(fc)], the Fourier space 'wiggles only' and the 
spherical harmonic correlation function [C(^)]. 

In section[2]we overview the different features which are 
present in the galaxy power spectrum, and discuss the po- 
tential information carried by each feature. In section [3] we 
describe the three BAO Methods used in this paper. In sec- 
tion |4] we describe the details of the Fisher forecast method 
used, as well as assumptions about our fiducial cosmologi- 
cal model and the implementations of each method. In this 
section we also describe the ground and space based sur- 
veys corresponding to Stages HI and IV for the BAO sur - 
veys described by the DETF report l|Albrecht et al.ll2006l ). 
In section [5] we compare our results with the DETF imple- 
mentation. We then present the constraints derived for each 
method for both fiducial BAO surveys and we also com- 
bine these with CMB constraints derived from the future 



^ http:/ /sci.osa.int/science-e/www/area/index.cfm?fareaid=102 



Planck mission. In addition, we investigate the impact of 
the wavenumber range on the Figure of Merit for each BAO 
method. In section[6]we discuss the constraints obtained by 
each method, and the hierarchy that exists between them 
and give our overall conclusions. Details of the implementa- 
tion of the Fisher matrix calculations for each BAO method 
as well as the for the Planck mission are given in Appendices 
[Handini 



2 THE BUILDING BLOCKS OF THE GALAXY 
POWER SPECTRUM 

The Fourier space matter power spectrum describes the fluc- 
tuation of the matter distribution and is defined by: 



(1) 



where 5{k) represents the Fourier transform of the matter 
overdensities S{f) — , and the mean density of the 

Universe is p. The term S'^{k~k') represents the Dirac delta 
function. 

The galaxy power spectrum is a rich statistic, where 
several features on different scales contain specific cosmo- 
logical information. On linear scales, we identify three main 
features in the galaxy power spectrum, which are: 

• The broad-band power: Information is contained in the 
shape, normalization and time evolution of the power spec- 
trum. 

• The baryon acoustic oscillations: Information is con- 
tained in the tangential and radial wavelengths, as well as 
the wiggle amplitude. 

• The linear redshift space distortions 

The three building blocks of the observed galaxy power 
spectrum are plotted separately in Figure [T] (panels B, C 
and D), while the observed power spectrum is plotted in 
panel A. Each building block probes the various sectors of 
the cosmological model (dark energy, dark matter, initial 
conditions) in a different way. In this section, we describe 
the main building blocks of the power spectrum and focus on 
their cosmological parameter dependence for the following 7 
parameters: [Qm, ^t, f^DE, ^fo, Wa, h, ris]. 



2.1 Broad Band Power 

The initial dark matter power spectrum is assumed to be 
of the form: Pinit(fc) oc k"" . The spectral index Us controls 
the tilt of the initial power spectrum, and the Harrison- 
Zel'dovich value of = 1 corresponds to a scale invariant 
power spectrum of perturbations. The linear matter power 
spectrum at a redshift z can be obtained by assuming the 
initial power spectrum has evolved according to: 



P^{k,z)=D\z)Pinit{k)T^{k). 



(2) 



The quantity T{k) is the transfer function and is a solution 
to Boltzman's equations which includes the baryon wiggles. 
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Figure 1. Illustration of the main building blocks of the linear galaxy power spectrum. Panel A: The observed linear galaxy power 
spectrum. This includes the linear galaxy bias, the smooth power spectrum (broad band power), baryonic wiggles and linear redshift 
distortions. In linear theory the redshift evolution of the galaxy power spectrum depends solely on the linear growth factor (see equation 
[2]|; we illustrate this by comparing the linear galaxy power spectrum at redshift z = (solid blue line) and redshift z = 0.5 (dot-dashed 
red line). Panels B, C and D illustrate the different building blocks of the galaxy power spectrum. Panel B: The smooth part of 
the galaxy power spectrum contains information in the shape and normalization of the power spectrum; it can be calculated using the 
analytic formula given by Eisenstein & Hu (1999). Panel C: The ratio (blue solid line) of the full power spectrum and the smooth part 
of the power spectrum reveals the residual baryonic wiggles. The dashed line corresponds to no baryonic wiggles. Panel D: The ratio 
(blue solid line) of the radial galaxy power spectrum (/i = 1 in equation to the tangential spectrum (/i = 0) illustrates the scale- 
independent effect of linear redshift distortions. Linear redshift distortions add power in the radial direction. The dashed line corresponds 
to no redshift distortions. 



The linear growth factor D{z), normalized to 1 today, quan- 
tifies the evolution of linear growth of structure with red- 
shift. Both the matter Qrn and the dark energy floE densi- 
ties affect the normalization of the power spectrum. 

The linear power spectrum today is characterised by a 
turnover (which occurs around k^q ~ 0.015/i~^Mpc in Fig- 
ure [!}; the turnover corresponds to the size of modes which 
enter ed the horizon at matter-radiation equaUty (jPodelso n 
l2003l ). The scale feoq at which this occurs is proportional to 
ilmh^ (with no dependence on the radiation density, since 
this is well constrained by the measurements of the CMB 
temperature anisotropies) . 

In a flat universe, the dark energy density f^DE will in- 
directly affect the position of the turnover through Ode ~ 
1 — ilm- In this work we always consider a curved universe 
where Qde and Sim are independent of each other, so that 
the dark energy density does not affect the position of the 
turnover but it does affect the overall amplitude. In the case 
where dark energy clusters (i.e. its equation of state param- 
eters are given hy wo 7^ —1 and Wa ^ when using the 
parameterisation defined in equation I13[l. the power spec- 
trum will be enhanced on large scales I Ma et al]|l999l ). 

Increasing the baryon density, Oj,, will have the efi^ect 
of decreasing power on small scales as well as also introduc- 



ing baryon wiggles in the power spectrum (see section . 
In this paper, the matter density is the sum of the Cold 
Dark Matter (CDM) and baryon (b) densities respectively, 
i.e. Qrn = f2cDM + Qb, where we choose Qm and to be 
independent. This means that when increasing Qb, Qcdm 
must decrease in order for f2,„ to be kept fixed. 

In a galaxy survey, the observable is the galaxy over- 
density (5g which is assumed to trace the underlying matter 
distribution 5m following: 



5g = 6(z, k)5„ 



(3) 



where b{z, k) is the galaxy bias, which can be both redshift 
and scale dependent. The real space galaxy power spectrum 
is then given by 



Pg(k,z) = b\z,k)D'\z)Pinitik)T\k). 



(4) 



On linear scales, it is often assumed that the scale 
dependence of the galaxy bias can be dropped so that 
b{k,z) — b(z); in this case, the linear bias only modulates 
the overall amplitude of the galaxy power spectrum Pg{k). 
The bias at redshift zero [b{z — 0)] is degenerate with the 
normalization as, so that in general the bias b{z) is degen- 
erate with both the normalization and the growth ■ D{z). 
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2.2 Baryon Acoustic Oscillations 

The baryon acoustic oscillations in the CDM power spec- 
trum are relatively weak features (see Figure never- 
theless t hey can be us ed independe ntly to c onstrain cos- 
mology (|Peacock et al.ii 2006; Albrc cht et alj [2006 ). They 
can be isolated from the galaxy power spectrum either by 
taking the ratio with the correspondi ng baryon-free power 
spectrum, or a smooth fittin g curve (|Blake fc Glazebrookl 
l2003l : ISeo fc EisensteinI [20031 ) . In the latter case, the oscil- 
lations can be approximated by a decaying sinusoidal func- 
tion. The amplitude of the oscillations increases with the 
baryon density Q.^,, and the location kA of the peaks is re- 
lated to the sound horizon at decoupling, s, by kA = 27r/s 
iBlake fc Glazebrookl [2OO3I ). 

The theoretical value of kA provides a known - or stan- 
dard - ruler, fixed by the sound horizon at decoupling: 



—At = — 

a Hq 



E{z) 



-dz, 



(5) 



where Cs is the sound speed and: 



E{z) = H{z)/Ho = v/aad(l + z)^ + + 2)3 + nuE.f{z) 



where the expression 



exp 



3/ i+^dz' 
1 + z' 



(7) 



describes the effect of dark energy on the Universe's expan- 
sion. At these high redshifts the curvature can be neglected. 
The radiation density is strongly constrained by the mea- 
surement of the CMB temperature, and the sound horizon 
s depends strongly on Qm through equation [6l In addition, 
the sound speed and the redshift of decoupling ^dcc depend 
on flbh^ , which is thereby moderately constrained, but apart 
from this combination, the constraints on Qb and Ho sepa- 
rately are weak due to degeneracies (though including the 
amplitude of the baryon wiggles would further constrain Qb)- 
Dependence on dark energy parameters can vary with 
the model, but in most cases the dark energy density at 
very e arly times is small and can be neglected ([Bean et al.l 
l|200ll ) showed it has to be less than a few percent during 
the Big Bang Nucleosynthesis epoch). The dependence on 
the dark energy parameters does therefore not enter through 
the size of the ruler, it affects the observed oscillation scale: 
A ruler at a redshift z with a given (comoving) size rj_{z) 
perpendicular to the line of sight, and a size r^^{z) along the 
line of sight is related to the observed angular size A9 and 
redshift extent Az through (Se o fc Eisenstein.,2003i ') : 



r±iz) = 



{l + z)DA(z)Ae, 
cAz 



H(z) 



(8) 

(9) 



where Da{z) is the angular diameter distance. Since the 
size of the ruler scales with s, the observation of the baryon 
acoustic scale in the tangential and radial direction affords 
us a measurement of Da/s and H-s respectively. In the limit 
where the size of the sound horizon is known from CMB 
data to much higher precision than the measurement of the 
oscillation scales from the galaxy survey, we can consider the 



baryon wiggles to provide estimates of the angular diameter 
distance and the expansion rate directly. 

Thus, if we divide out the broad shape of the power 
spectrum and neglect the amplitude of the wiggles, we ex- 
pect to be able to measure the matter and dark energy den- 
sities Qrn and Qde', with sufflcient amount of redshift infor- 
mation, it will also be possible to constrain the evolution of 
the dark energy equation of state w{z). 

2.3 Redshift Space Distortion 

An observer can only measure the galaxy power spectrum 
in redshift space, which is distorted compared to the power 
spectrum in real space. Redshift distortions are due to pe- 
culiar velocities of galaxies; these cause radial distortions 
in the observed galaxy density field. This distortion occurs 
because the observed redshift Zobs is a sum of two quantities: 



Zobs 



Zh ~t- 2t7. r, 



(10) 



where Zh is the redshift due to the cosmological Hubble ex- 
pansion of the Universe, and z^,p is the redshift due to the 
j-g-j radial component of the galaxy's peculiar velocity. 

In linear theory this will affect structures along the line 
of sight, which will appear enhanced compared to transverse 
structures; i.e. for a structure which is isotropic in real space, 
an observer will measure more power in the radial direction 
than in the transverse direction. Considering linear theory 
alone, the observed redshift space power spectrum for galax- 
ies is related t o the real space matter power spectrum by 
(Kaiser 1987; Seo fc Eise nstein 2003): 



'\k,z,^) oc (l-f /3^2)%(z,fc)2 



\k,z) 



(11) 



where ^ is the cosine of the angle between the wavevector k 
and the line of sight. The amplitude of the redshift distortion 
is modulated by the distortion parameter 



1 dlnL»(a) 
b{z) dlna 



(12) 



so that the redshift distortions are a probe of the growth 
rate of structure as well as the linear galaxy bias. 



3 RECIPES FOR BARYON ACOUSTIC 
OSCILLATIONS 

The goal of this paper is to contrast the constraining poten- 
tial of three different families of BAG methods, namely: 

(i) Full Fourier space galaxy correlation function P{k), 

(ii) Fourier space BAG 'wiggles only', 

(iii) Spherical harmonic space galaxy correlation function 
C(^). 

The general underlying basis of each method is de- 
scribed in this section and a schematic summary is given 
in Table [1] Details of the tomographic method, the fiducial 
cosmological model and further implementation specifica- 
tions are given in section [l] Full mathematical details on 
the Fisher matrix calculations for each family of methods 
are given in Appendix 1X1 
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Table 1. Schematic summary of which building blocks of the galaxy power spectrum are probed by each family of BAO method 
considered in this paper. In this paper, certain components such as the normalization of the power spectrum or the redshift distortions 
have been marginalized over in order to follow prescriptions which exist in the literature; when this is done the corresponding building 
block is represented in brackets in the Table below. 





Broad Band 
Power 


Overall Amplitude 


Tangential 
BAO scale 


Radial 
BAO scale 


Wiggle Amplitude 


Redshift Space 
Distortions 


Full P{k) 




W) 








W) 


'Wiggles only' 


















W) 











3.1 Full P(k) 

The full Fourier space galaxy correlation method uses in- 
formation across different scales and capitalizes on each 
building block of the galaxy power spectrum. Such meth- 
ods have already been applied to forecasts for future surveys 
l|Seo &: Eisensteinl2003l : lAmendola et al.ll2005l . see Appendix 

In this paper w e follow the implementation of 
ISeo fc EisensteinI l|2003l l. which we refer the reader to for 
specific methodology, though details are also given in sec- 
tion 2] and Appendix I All 

Intuitively, because this method uses information from 
the full galaxy correlation function, it should have the po- 
tential to constrain cosmological parameters with high pre- 
cision. However since it uses information over a wide range 
of scales this method could also be prone to a high level of 
systematics. 

In particular, the unknown linear bias b{z) will affect 
both the overall amplitude of the power spectrum as well as 
the amplitude of the redshift distortion (see equation [11]). If 
the bias is also scale dependent, i.e. if b{z) — h(k,z), then 
this can potentially distort the shape of the power spectrum 
over a range of scales leading to a fairly high sensitivity to 
systematics. 

The full galaxy correlation function will also be more 
sensitive to non-linear redshift distortions at small scales. 
Additionally, intrinsic non-linearities in the density contrast 
will distort the shape of the spectrum at large wavenum- 
bers - while leaving the w i ggle lo catio n relatively unchanged 

i see fo r e.g. IScoccimarrol 1 2004 ) and iMatarrese fc Pietronil 
200i)). 



3.2 BAO wiggles only 

To avoid the effect of the potential systematics of the full 
P(fe) method, one can focus the analysis on specific scales, 
and remove the overall shape and amplitude of the galaxy 
power spectrum, centering the analysis on the baryon wig- 
gles only. 

The way to do this is to consider the ratio of the 
observed power spectrum to a baryon-free power spec- 
trum, i.e. P{k,Q,b 7^ 0)/P{k,Qi = 0) (as in Figure 
1]) or to a smooth p arametric curve l|Blake et al.l l2006l : 



Seo fc Eisenstein|[2007l ). This BAO 'wiggle-only' method has 



been implemented in several papers. Some popular methods, 



are described inlParkiiison et all (|2007l ') , Isiake et all (|2006l ) 
and lSeo fc EisensteinI l|2007h . 

By focusing on the baryon wiggles only, the measured 
quantity is now independent of the redshift dependent lin- 
ear bias b{z) (though the errors on the peak measurement 
will depend weakly on the bias) and, because the wiggles oc- 
cur on a limited wavenumber range, the measured quantity 
is also weakly dependent on the scale dependent bias b{k). 
This method should intuitively be more robust than the full 
P{k) approach, as it exploits less features over a limited k 
range, it is also bound to provide weaker constraints on the 
cosmological parameters. 

In this paper, the main calculations for the 'wig- 
gles only' metho d are perfor med using the f o rmali sm of 
ISeo fc EisensteinI l|2007l ') and iParkinson et al.l (|2007l ) (see 
section |3] and Appendix IA2[) . except in section [5.11 where 
the errors on the wiggles estimation are calculated as in 
iBlake et al.l l|2006l ). 



3.3 Spherical Harmonics C{£) 

Both methods described above perform the analysis in 
Fourier space. Doing this raises several issues. 

The first is that equation (|ll|l which relates the real 
and redshift space Fourier power spectra is only valid in 
the far field approximation, i.e. for galaxies which are sep- 
arated by small angles on the sky, and for which the line 
of sight vectors can be considered parallel. While this may 
be a valid approximation for current surveys, future sur- 
veys plan to cover all the extragalactic sky (~ 27r) and 
this may no longer be valid. A natural decomposition for 
data on a sphere is in spherical harmonic space in which 
transverse and radial modes are independent. In this de- 
composition, there exists an exact solution to ex press the 
galaxy power spectrum in red shift space (see iFisher et al.l 
Il994l : [Heavens fc TavloJ[l995l . and equations IXTOl and [XTOl 
in Appendix 1^. For large sky coverage it is therefore more 
judicious to study the galaxy power spectrum in spherical 
harmonic space. 

The second and more fundamental issue has to do with 
the actual measurement of the Fourier power spectrum. The 
bare observables provided by a redshift survey are the angu- 
lar position and redshift of a galaxy, say {6, 0, z) (although 
technically the redshift is itself a first order quantity derived 
from the galaxy's observed electromagnetic spectrum). To 
measure the Fourier space power spectrum it is necessary 
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to relate these observables to the wavenumber k and this 
transformation is dependent on the assumed cosmological 
model. On the other hand, in spherical harmonic space, the 
transform from (6, (j), z) to (I, z) is independent of cosmology 
and this is a further central motivation for analysing galaxy 
surveys in spherical harmonic space. 

This is the final family of methods we in vestigate, fol- 
lowing the tomographic method described in iDolnev et al.l 
(t2006 ) (see section |4] and Appendix I A3I for details). 



4 IMPLEMENTATION DETAILS 

In this section we describe the implementation details of our 
calculations. In section \AA[ we justify our choice of fiducial 
cosmological parameters and their central values. In section 
14.21 we describe the two fiducial surveys we use for our cal- 
culations, which are chosen as Stage III and Stage IV BAO 
surve ys as described in the DETF report ijAlbrecht et al.l 
I2OO6I ): these correspond to future ground and space based 
surveys respectively. In section 1431 we overview the specifi- 
cations of our calculations. 



4.1 Fiducial cosmological model and central 
values 

To quantify the dark energy, we adopt a convenient param- 
eterisation ijChevaUier fc Polarskill200ll : lLinder..2003 ) of its 
equation of state w{a) by. 



10(a) 



(13) 



where a = l/(l + z). This can also be expressed in terms of 
the pivot scale factor a^: 



10(a) = Wp + {ap — a)iOa. 

The DETF "Figure of Merit" (FoM) is defined by: 
1 



FoM = 



a{wp)a{wa) ' 



(14) 



(15) 



which is used to compare the suitability of different probes 
and surveys to constrain the dark energy equation of state. 
In order to avoid problems of phantom crossing around wq = 
— 1, we chose a central value of wq — —0.95. 

We forecast errors for an 7-parameter model within the 
framework of general relativity, i.e. we have not added any 
extra para meters to account for potential modified gravit y 
scenarios l|Amendola et all I2OO8I : iHuterer fc Linde j l200'i1 ) . 
The 7-parameter model and central fiducial values are given 
in Table [21 The power spectrum normalization is chosen as 
(jg — 0.80 and is marginalised over for calculations involving 
the full Fourier power spectrum and the spherical harmonic 
galaxy correlation function. 



4.2 Fiducial Survey 

The DETF have published dark energy forecasts for differ- 
ent "Stages" of surveys, numbered from I to IV. Stage I 
corresponds to dark energy projects that have already been 



Table 2. Central fiducial values and cosmological parameter set 
for the Fisher matrix calculations performed in section [5] The 
normalization of the power spectrum is taken as a-8=0.80 . 



Parameter Central fiducial value 



Wo 

h 

ris 



0.25 
0.0445 
0.75 
-0.95 


0.70 

1 




2.0 



Figure 2. The value of fcmax as a function of redshift z. The 
solid line (diamonds) corresponds (j{R) < 0.20 and fcmax < 
0.25/iMpc~^. The dashed and dot-dashed lines correspond respec- 
tively to deviations around the central line, specifically fcmax (^) • 
0.5 (dot-dashed) and fcmax (2:) ■ 1-2 (dashed) respectively and are 
used in section [5.31 



completed; Stage II to ongoing projects. Stage III to "near- 
term, medium-cost, currently proposed projects" and Stage 
IV to future ground or space based (LSST/SKA or JDEM 
like) missions with large sky coverage. The BAO forecasts 
published by the DETF were calculated using a specific im- 
plementation of the 'wiggles only' method. In this paper we 
are interested in comparing these forecasts with those from 
two other methods, namely the full Fourier power spectrum 
and the spherical harmonic power spectrum. 

We make these comparisons for future ground and space 
based spectroscopic ("S" identifier in DETF denomination) 
surveys, which correspond to "Stage III" and "Stages IV" 
surveys. In the DETF report, the spectroscopic Stage III 
survey is composed of both a wide and a deep survey 
(BAO-IIIS-o in DETF denomination, defined page 56 in 
lAlbrecht et aP (|2006l )). for which results are always pre- 
sented combined. Here we consider only the wide survey of 
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Table 3. Future fiducial spectroscopic ground and space based surveys taken from DETF. The sky coverage, redshift range and galaxy 
density (nP = 3) are described in the DETF report. The number of redshift bins, the bias prescription and wavenumber range are chosen 
specifically for this paper. 





Ground Based Survey 


Space Based Survey 


DETF denomination 


BAO-IIIS-o (WFMOS-like, "wide" only) 


BAO-IVS-o (Optical/NIR JDEM Spatial Mission) 


Sky coverage 


2, 000dcg2 = 0.05) 


10,000deg2{/3ky = 0.25) 


Redshift range 


0.5 < 2 < 1.3 (4 bins) 


0.5 < 2 < 2 (8 bins) 


Bias prescription 


feg(^) = Vi 


+ z 


fe-rangc 


fcmin = 10-=^, fcmax < 0.25 Or Ct{R) < 0.20 



the Stage III. For the Stage IV mission, we chose it to resem- 
ble a spectroscopic JDEM-like mission, by u sing the DETF 
notati on BAO-IVS-o (defined page 56 in lAlbrecht et al.l 
l|2006t ) - the optimistic survey does not include systemat- 
ics). 

The main features of these two spectroscopic surveys 
are given in Table O The galaxy number density was taken 
so that neffPg(fc* = 0.2/iMpc~^) = 3 (| Blake et al.l l2006l : 
lAlbrecht et "al]|2006h . which fixes the galaxy distribution and 
shot noise contribution. We note that fixing the value of 
nesP{k*) does not construct a realistic magnitude-limited 
galaxy distribution. The sky coverage and redshift ranges 
are chosen according to the description given in the DETF 
report and given in Table [S] Specifications given in the sec- 
ond half of the Table are chosen specifically for this paper 
and are explained in more detail below. 



4.3 Specifications 

For the 'wiggles only' calculation (which is the calculation 
performed by the DETF) the value of the galaxy bias is ir- 
relevant - so long as it is scale independent, which is what we 
assume here. For the Fourier space and spherical harmonic 
space calculations, we chose the bias to follow a simple ad 
hoc functional form of b{z) = \/l + z. This choice can po- 
tentially affect statistical constraints for the full P{k) and 
C{£) methods, where the bias is explicitly present; whereas 
the 'wiggles only' method depends very weakly on the value 
of the galaxy bias. 

As the linear bias and normalization erg are highly de- 
generate, we also marginalise over the bias value b{zi) in 
each redshift bin as well as over as for the C(i) me thod 
(which is slightly different than the lDolney et aL |(|2006) im- 
plementation, who uses a more complex bias prescription). 

For the full Pjk) met hod we follow the implementation 
of ISeo fc EisensteinI (|2003h : information from the amplitude 
of the power spectrum as well as from redshift distortions are 
suppressed by marginalizing over the product asb{zi)D(zi) 
and the redshift distortion parameter /9(zi) in each redshift 
bin. Finally, we margina lize over an unknown w hite shot 
noise, which is denoted in lSeo fc EisensteinI (|2003l ) by Pshot- 

For all methods we perform calculations in the linear 
regime. This is ensured by restricting the calculation to cer- 
tain linear scales. The fcmax at which non-linear effects be- 
come no longer negligible in our fiducial model depends of 



course on the precision one wishes to attain. Recen t work 
(jCrocce fc Scoccimarrd l2008l : iMatarrese fc Pietronil |2008| ) 
show that at A: = 0.2/iMpc~^ the deviation from linearity at 
z—0 is already sizeable, although the position of the peaks 
is much less sensitive to non-linear corrections. 

The wavenumber corresponding to the linear cut-off will 
evolve with redshift. We quantify this evolution by con- 
sidering only scales for which (t(_R) < 0.20 and fcmax < 
0.25/!,Mpc~^, where o"(-R) is defined similarly to the nor- 
malization as = a{R = 8/i~^Mpc), but for a general R. 
The evolution of the fcmax we use for our computation with 
redshift bin is illustrated in Figure [2] as the solid line with 
diamonds. We assume the largest scales probed are given by 
k^i^ = 10"^/iMpc"^ 

We find the choice of wavenumber range can have a very 
strong influence on results. We investigate the dependence 
of the FoM on the values of fcmax in section 15.31 



5 FORECASTS 

The 'wiggles only' BAO method we use is comparable to 
that used by the Dark Energy Task Force (jAlbrecht et al.l 
|2006| . DETF), and so we first compare constraints from our 
calculation with those of DETFas10 in section 15.11 In sec- 
tion [S3] we compare the constraining potential of the three 
different families of BAO methods. The constraints are cal- 
culated for the galaxy surveys described in section|4l namely 
a future ground- and space-based spectroscopic surveys. 

5.1 Comparison with DETF Calculation 

In this paper, constraints for the BAO 'wiggles only' 
method are performed using distance error estimates from 
ISeo fc EisensteinI (|2007l ). The DETF uses a similar method, 
though the distance error esti mates are ca l culate d us- 
ing the fitting formulae given by iBlake et al.l (|2006l ). The 
slight differ e nce b etween the fitting formulae provided in 
Blake et al.l bOOdt) is th at the fc-range is fixed, whereas in 
Seo fc EisensteinI l|2007l ) the fe-range is an input to the er- 
ror estimates. In this section we sw itch to using the er- 
ror estimates from iBlake et al.l (|2006l ) in order to compare 

^ http://www.physics.ucdavis.edu/DETFast/ by Jason Dick and 
Lloyd Knox. 
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constraints with those given by the DETF. We obtain con- 
straints from DETFast^ , a web-interface for calculations per- 
formed by the DETF. 

In this calculation the central fiducial model and cosmo- 
logical parameter set was changed slightly to adopt the same 
as that in the DETF report, i.e. Q,de = 0.7222, Q,k = 0, 
WQ — —1, Wa = 0, and Qmh^ = 0.146 and the sound horizon 
is kept fixed. The fiducial survey was chosen to be BAO-IVS- 
o, as described in Table |31 Comparison between constraints 
is given in Table O the reported constraints are shown with- 
out any priors (i.e. without Planck priors). Both calculations 
show agreement at the sub-percent level for all marginalised 
errors (as well as for fixed errors - not shown in this table). 



5.2 Comparison of different BAO Methods 

For the comparison between the three different BAO meth- 
ods, we switch back to the fiducial cosmological model 
given in Table |2[_ We also use error estimates from 



ISeo fc EisensteinI (|2007l ) for the 'wiggles only' method, in 
order to have control over the fc-range. 

The three methods are: the full power spectrum in 
Fourier space [P{k)], the Fourier space wiggles only method 
[BAO 'wiggles only'] and the full projected power spectrum 
in spherical harmonic space [C(^)]. In TableUwe present the 
marginalised cosmological constraints for the 7-parameter 
model we consider from these three different BAO methods. 
These constraints are marginalised over the linear bias and 
power spectrum normalization (as) for the C{£) method; 
for the full P{k) method we marginalise over the amplitude 
of the power spectrum and the redshift distortions term as 
detailed in section [473] 

The results in brackets have been combined with fore- 
casted constraints from the Planck CMB mission, or what 
we refer to from now on as Planck priors. The forecasted 
CMB priors include scalar temperature and E-mode polar- 
ization information. Details on how we calculated the Planck 
priors are given in Appendix iBl 

Parameters for which the table entry shows ' — ', are not 
constrained by the corresponding method (i.e. the affected 
parameter is omitted from the Fisher matrix calculation), 
and we have not quantified constraints which were weaker 
than 10. The cosmological constraints for different methods 
are presented in Figure [3] for the dark energy parameters wo 
and Wa for the ground (top panel) and space (lower panel), 
without Planck priors. 

We first explore the results from Table 3] without Planck 
priors as this gives us an understanding of what each BAO 
method constrains. 

We begin by exploring the ground-based survey with- 
out Planck priors. The 'wiggles only' method, which uses no 
information from the broad band power, the redshift space 
distortions nor the wiggle amplitude (see Table[T|, gives con- 
straints which are up to a factor of 7 weaker than the C{£) 
method. Furthermore, the 'wiggles only' method fails to con- 
strain either Q^, Wa, h and Us, whereas the C{£) method only 
fails to constrain Wa- The constraints from the C{£) method 
are a factor 2 to 18 weaker than those from the P{k) method, 
which can constrain Wa (though the errors are still large). In 



all cases, the Figure of Merits (FoMs) are small, making any 
comparison between them meaningless (for very low FoMs, 
upper limits only are shown in Table |4}. 

Adding the Planck priors to the ground-based case im- 
proves constraints (hereafter we omit improvements on pre- 
viously unconstrained parameters) from the 'wiggles only' 
method by a factor of 3 — 13, from the C{£) method by a fac- 
tor of 1.5— 10 (omitting the improvement by a factor of 300 
on Us), from the P{k) method by a factor of 2 — 37. Includ- 
ing the Planck priors also reduces the differences between 
the 'wiggles only' and the C{£) methods (the marginalised 
constraints are of the same order of magnitude for all pa- 
rameters), but the constraints from the P{k) method are for 
some parameters up to a factor 16 tighter than for the C{£) 
method. The FoM for the P{k) method is over a factor of 
17 larger than the FoM for both the 'wiggles only' and C{£) 
method. 

We now consider constraints for the space-based survey, 
beginning with the constraints without priors. 

For the space-based survey, constraints from the 'wig- 
gles only' method are a factor of 1.5 — 5 weaker than those for 
the C{£) method, and fij, and h still remain unconstrained 
for the 'wiggles only' method (which is expected, since this 
method constrains only the product Qmh^, see l2.2|l . 

The constraints from the P{k) method remain consid- 
erably stronger, giving an improvement of 2 — 21 on both 
the 'wiggles only' and the C{£) method, where the largest 
improvement (factor of improvement greater than 15) is on 
h, Us, 0.m and wo. The FoM for the P{k) method is 52, i.e., 
35 times larger than that for the 'wiggles only' method and 
12 times larger than that for the C{t) method. 

Adding the Planck priors to the space-based case im- 
proves constraints from the 'wiggles only' method by a factor 
of 3 — 5, on the C[t) method by a factor of 1.5 — 7 (omit- 
ting the factor of 40 improvement on n^), and on the P{k) 
method by a factor of 1.3 — 14. As with the ground-based 
survey, the space-based survey with the Planck prior shows 
little variation between the 'wiggles only' and the C{£) meth- 
ods, with the Figure of Merit for each method being the same 
order of magnitude. The FoM for the P(k) with Planck pri- 
ors is 502, approximately 18 times larger than the FoM for 
the 'wiggles only' method with Planck priors, and a factor 
of 30 tighter than C(^)' method with Planck priors. 

The constraints for wo, Wa, f^DB and Wp for the space- 
based survey, using the 'wiggles only' method with Planck 
priors are of the same order of magnitude of those pre- 
sented by the DETF for BAO-IVS-o ~f Planck (page 77 of 
lAlbrecht et al] (|2006l )). though comparison should be taken 
with caution as both use different calculations for the Planck 
prior as well as different methods of calculating the errors 
for the BAO calculation. The DETF report does not quote 
constraints on the other parameters, though they have also 
been marginalised over. 

In order to investigate which building block of the 
power spectrum provides the most information to the P{k) 
method, we evaluate new constraints, this time for a quasi 
baryon-free spectrum (to do this we set ^Ib = 0.005 as lower 
values causes numerical problems as other authors have also 
found). This test suppresses the baryon wiggles, leaving the 
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Table 4. Marginalised Icr constraints from the three different BAO Methods for DETF like ground-based survey (BAO-IIIS-o) and 
space- based survey (BAO-IVS-o). 





Fourier 


space BAO 'wiggles only' 


Spherical Harmonic Space C{t) 


Fourier 


space Full P{k) 


Ground Based Survey (with Planck priors) 




2.04 


(0.170) 


0.303 


(0.162) 


0.137 


(0.01) 




> 10 


(0.0303) 


0.122 


(0.0290) 


0.0310 


(0.002) 




1.99 


(0.156) 


1.25 


(0.163) 


0.215 


(0.016) 


WO 


6.28 


(2.42) 


3.23 


(1.69) 


0.56 


(0.25) 


V]p 


1.15 


(0.184) 


0.979 


(0.650) 


0.358 


(0.07) 


Wa 


> 10 


(5.41) 


> 10 


(4.45) 


3.07 


(0.65) 


h 


> 10 


(0.238) 


2.38 


(0.227) 


0.13 


(0.016) 






(0.00465) 


1.54 


(0.00465) 


0.165 


(0.0045) 


Dark energy FoM 


< 0.090 


(1.2) 


< 0.10 


(0.35) 


0.91 


(21) 


Space Based Survey (with Planck priors) 




0.181 


(0.0385) 


0.0367 


(0.0254) 


0.00985 


(0.003) 




> 10 


(0.00687) 


0.0193 


(0.00463) 


0.00221 


(0.0006) 




0.178 


(0.0375) 


0.113 


(0.0251) 


0.0692 


(0.005) 


WO 


0.990 


(0.377) 


0.355 


(0.331) 


0.0638 


(0.049) 


Wp 


0.197 


(0.0410) 


0.147 


(0.0955) 


0.0313 


(0.02) 


Wa 


3.94 


(0.878) 


1.63 


(0.610) 


0.612 


(0.10) 


h 


> 10 


(0.0539) 


0.244 


(0.0360) 


0.0118 


(0.004) 


rig 




(0.00463) 


0.182 


(0.00450) 


0.0106 


(0.0034) 


Dark energy FoM 


1.5 


(28) 


4.2 


(17) 


52 


(502) 



Table 5. Comparison of Icr marginalised constraints from BAO 
'wiggles only' (without Planck priors with the DETFast results. 
The results agree at the sub-percent level. 



BAO 'wiggles only' DETFast 

0.0261 0.0258 

0.0971 0.0967 

Hfc 0.0308 0.0307 

wo 0.612 0.610 

Wp 0.120 0.119 

Wa 2.51 2.49 

Dark energy FoM 3.32 3.37 



broad band power as the only building block which now 
provides information, since the amplitude and the redshift 
distortions are marginalised over. 

For the space-based survey, we find a new FoM of 11 
(123) without (with) Planck priors, compared to the pre- 
vious values of 52 (502) for the calculation including the 
baryon wiggles. Though the values of the FoM for the P(k) 
method have dropped by about a factor of 4 when omitting 
information from the baryon wiggles, these constraints are 
still an order of magnitude tighter than those for the 'wiggles 
only' method, suggesting that a large part of the constrain- 
ing power of the P{k) method comes from the information 
contained in the broad band power. 




-6 -4 -2 2 4 

Wo 




-1.6 -1.4 -1.2 -1.0 -0.8 -0.6 -0.4 
Wo 



Figure 3. Marginalised Icr constraints (without Planck priors) on 
dark energy equation of state parameters wo and Wa for the three 
different BAO methods: Red (dotted): Full P(k), Blue (solid): 
C(l), Green (hashed: BAO 'wiggles only'. TOP: ground based 
survey (BAO-IIIS-o). BOTTOM: space-based survey (BAO-IVS- 
o). 



© 0000 RAS, MNRAS 000, 000-000 



10 Rassat et al. 



No Priors With Planci< Prior 





1 000.0 


1 1 1 1 1 1 1 1 1 1 1 1 1 


. . 1 1 


... 1 ... 1 ... 1 ... 1 1 
.,,.o 












g 


100.0 








< 










* 

a. 






















10.0 




























FoM 


1.0 


/ y 




P(k) 






/ y 
/ y 
□ 

y 




-B- C(l) 






y 

y 

a' 




Wiggles Only 




0.1 


... 1 ... 1 ... 1 . 


. . 1 , 


... 1 ... 1 ... 1 ... 1 , 



0.4 0.6 O.B 1.0 1.2 0.4 0.6 O.B 1.0 1.2 

Normolized k..^, Normolized 



Figure 4. Figure of Merit (FoM) as a function of normalized fcmax 
for the full power spectrum P(fe) [dotted, diamonds], the spherical 
harmonic C(€) [solid, squares] and the 'wiggles only' [dot-dashed, 
triangles] methods - for the space based survey BAO-IVS-o. The 
normalized values of fcmax a.re taken from Figure [2] The left-hand 
panel corresponds to constraints from the BAG methods without 
any priors, the right-hand panel includes CMB constraints from 
Planck. 



5.3 Dependence of Figure of Merit on 
wavenumber range 

We find the FoM can depend strongly on the choice of non- 
linear cut-off or fc-range. We investigate this dependence by 
considering three different maximum wavenumbers (keeping 
fcmin fixed at 10~^/iMpc~^), namely fcmax(2)-0. 5, k-^^-^{z)-\.2 
and the previous value of fcmax(2). These three wavenumber 
ranges evolve with redshift and have been plotted in Figure 

m 

The evolution of the FoM is shown for the three meth- 
ods considered (full P{k), spherical harmonics C{i) and the 
'wiggles only') for the space-based survey (BAO-IVS-o) in 
Figures |4] and [5] The left hand panels in both show the re- 
sults without any priors and the right hand panels include 
Planck priors. The FoM for the P{k) method evolves rapidly 
with fcinax, jumping from 3.2 at fcmax(«) • 0.5 to 51 at k^^^{z) 
and 80 at fcinax(2) ■ 1-2. The evolution of the C{£) method is 
smoother: from 0.4 at fcmax(z) • 0.5, to 4.2 at k^^^^(z), and 
7.1 at fcmax(z) ■ 1-2. Finally for the 'wiggles only' method, 
the evolution is slowest: from 0.2 at fcmax(z) ■ 0.5, to 1.7 at 
k^axiz) and 2.4 at fcniax(«) ■ 1-2. In Figure O we normalize 
the FoM to its value at km,ax{z) ■ 0.5 in order to quantify 
the factor of improvement of the FoM as a function of fcmax- 
This is plotted in Figure [S] with the left panel corresponding 



to results without any priors, and the right panel to results 
including the Planck CMB priors. 

Using the normalized values, we clearly see the rela- 
tive trends for each BAO method as a function of kmaxiz). 
Whether we consider the results with no prior or those in- 
cluding the Planck CMB priors, both the P{k) and C(^) 
methods evolve more rapidly than the 'wiggles only' method, 
and the difference is especially noticeable as we include 
larger values of fcmax- 

The full power spectrum P{k) method contains infor- 
mation about the broad band power spectrum, and the 
amount of information used from the broad band power will 
depend on the maximum wavenumber included in the Fisher 
matrix calculation. 



6 DISCUSSION 

The BAO features of the galaxy power spectrum P{k) 
have recently been observed in the latest galaxy surv eys 
(|Eisenstein et al.ll2005l : lHiitsill2006al : iPercival et al]|07 al ibi). 
They are also predicted to become a fundamental tool 
for precision cosmol ogy in the future jPeacock et al.ll2006l : 
lAIbrecht et"al]|2006l ). However, there exists a range of BAO 
implementations in the literature. In this paper we discuss 
the relative information provided by the different methods, 
and quantify the potential of each method to constrain dark 
energy parameters. 

The three BAO methods we focus on are: 

• Full F ourier space galaxy correlation function P{k) (see 
ISeo fc Eisen stcin (2003)). 

• F o urier space BAO 'w iggles only' (see lSeo fc EisensteinI 
l|2007l ): lBlake et al.l (|2006l )). 

• Sphe rical harmoni c spac e galaxy correlation function 
C{£) fsee lDolnev et al.l (|2006l )). 

We do not consider the real space [^(r)] or configuration 
space [w{6)], correlation functions where, although the BAO 
signal is larger, the non- linear and linear scales are more 
difficult to separate. 

Each method probes a different combination of features 
in the galaxy distribution and these are summarised in Ta- 
bled The full Fourier space power spectrum method probes 
all of these features and is thus expected to provide tighter 
constraints than both the 'wiggles only' and C{£) method. 
The 'wiggles only' method probes both the radial and tan- 
gential BAO scales - where the C{£) method only includes 
information from the tangential scale (this is because the 
C{£) method uses the projected 2-dimensional power spec- 
trum). However the C{£) method includes information from 
the broad band power of the projected galaxy power spec- 
trum. 

We use the Fisher matrix formalism to forecast con- 
straints for two future fiducial surveys, which we choose as 
Stage III (ground-based) and Stage IV (space-based) sur- 
veys as described in (Albrccht et al. 2006, DETF). The spe- 
cific DETF denominations for these surveys are BAO-IIIS-o 
(wide survey only) and BAO-IVS-o. We compare constraints 
for a 7 parameter model which allows for spatial curvature. 
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The evolution with redshift of the dark energy equation of 
state w{z) is parameterised as in equation 1131 Our central 
fiducial model in given in Table[2l The survey and implemen- 
tation details are given in Table |3] The choice of non-linear 
cut-off is justified in section and its redshift dependence 
shown in Figure [21 

In Table Owe compare results from DETFas10 with our 
calculations and find agreement at the sub percent level, us- 
ing the same fiducial parameter set and values as the DETF. 

The constraints for each method are presented in Ta- 
ble [J] without any priors and we also combine these results 
with priors from the future CMB mission Planck (the Planck 
priors for the same cosmological model). 

The results in Table |4] show that there exists a hier- 
archy between the three different methods considered here. 
The P{k) method which uses information from each feature 
in the linear galaxy power spectrum gives the tightest con- 
straints for both the ground-based and space-based survey. 

The 'wiggles only' and C{£) methods use different in- 
formation in the galaxy power spectrum. The 'wiggles only' 
method uses both the radial and tangential modes, but does 
not use information from the broad band shape of the power 
spectrum, the redshift distortions, nor the amplitude of the 
baryon wiggles. On the other hand, the C{£) method does 
not use the radial scale of the BAO scale (this is because the 
C{£) method consists of the projected power spectrum). One 
could use the 3-dimensional decomposition of the galaxy 
field in spherical harmonic space in order to recover the ra- 
dial BAO information as well. 

Because of this, the hierarchy between the 'wiggles only' 
and the C{£) method cannot be predicted. Here we find 
that for the two surveys considered, the constraints for the 
'wiggles only' method are weaker than those from the C{£) 
method, themselves are weaker than constraints from the 
P{k) method. For the space based survey, the Figure of 
Merit (FoM, defined in equation llSp for the P{k) method 
is 35 times greater than that for the 'wiggles only' method 
and 12 times greater than for the C{£) method. 

We also tested how the FoM for the P{k) method 
changed when a baryon-free spectrum was used. The FoM 
dropped by a factor of 4, but was still an order of mag- 
nitude larger than the FoM for the 'wiggles only' method, 
suggesting that a large part of the constraining power of 
the P{k) method comes from information contained in the 
broad band spectrum. 

Adding the Planck CMB priors to those from each BAO 
method reduces differences in the hierarchy. In particular 
'wiggles only' constraints with Planck priors are the same 
order of magnitude as those for the C{£) method with Planck 
priors. In all cases constraints from the P{k) method are 
tighter (with or without CMB priors). 

We find the constraints can depend strongly on the cho- 
sen value of the non-linear cut-off, and so quantify this de- 
pendence for three choices of kmaxiz). The evolution of the 
FoM (with and without Planck priors) with fcmax(z) is plot- 



http:/ /www. physics. ucdavis.edu/DETFast/ by Jason Dick and 
Lloyd Knox. 
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Figure 5. Factor of improvement of the Figure of Merit (FoM) as 
a function of normalized fcmax for the full power spectrum P{k) 
[dotted, diamonds], the spherical harmonic C(€) [solid, squares] 
and the 'wiggles only' [dot-dashed, triangles] methods - for the 
space based survey BAO-IVS-o. The improvement is quantified 
relatively to the value of the FoM at fc^o™ = 0.5 ■ fc^^l^- The 
normalized values of fcmax are taken from Figure|2] The left-hand 
panel corresponds to constraints from the BAO methods without 
any priors, the right-hand panel includes CMB constraints from 
Planck. In the left panel (without Planck priors), it is clear that 
the FoM for the P(fc) methods evolves more rapidly with fcmax 
than the C(t) method, which in turn evolves more rapidly than 
the 'wiggles only' method. When Planck priors are included, the 
FoM from the C(t) method evolves more rapidly than the other 
two methods. This may be do to different parameter degeneracies. 
In this case, the FoM from the P{k) method evolves more rapidly 
than the that for the 'wiggles only' method at high fcmax. 



ted in Figure [4] We find that there also exists a hierarchy in 
the evolution of the FoM with the non- linear cut-off fcmax (z). 
In the case where we consider the FoM from the three differ- 
ent methods without external priors, we find that the FoM 
for the full power spectrum Pik) method depends strongly 
on fcmax, jumping from 3.2 to 80 for the extreme values of 
fcmax. The dependence of the FoM for the spherical har- 
monic method Cit) is weaker, from 0.4 to 7.1 over the same 
wavenumber range; finally, for the 'wiggles only' method, 
the dependence is the weakest: from 0.2 to 2.4 for the same 
wavenumber range. 

This suggests indeed that the 'wiggles only' method 
presents a more robust analysis, in that it depends weakly on 
the choice of non-linear cut-off and the wavenumber range. 
Both the -P(fc) and C(€) methods depend strongly on the 
wavenumber range, and will also depend on the linear bias 
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prescription. In any case, when discussing constraints from a 
BAO method, it is necessary to expUcitly state the method 
used, as we find the FoM can vary by up to a factor of 35 
between different methods. 

The main conclusions of this paper are therefore: 

• There is not a unique BAO method, but there exists a 
range of implementations. 

• For the three BAO implementations considered in this 
paper, we find that for a given fiducial survey and cosmology 
the figure of merit can vary by up to a factor of 35 between 
methods. 

• There exists a hierarchy in the constraining power of the 
methods. We find the more information used by a method, 

the higher its statistical performance, but the higher its sen- 
sitivity to systematics and implementations details. 
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APPENDIX A: 
FORECASTS 



DETAILS ON FISHER 



It is possible to forecast the precision with which a future 
experiment will be able to constrain cosmological parame- 
ters, by using the Fisher Infor mation Matrix (for a detailed 
deriv ati on of the follow ing see iTegmark. Tavlor fc Heavens! 
Il997l or lDodelsonll2003l ) . This method requires only three 
fundamental ingredients: 

• A set of cosmological parameters 8 for which one wants 



to forecast errors and a fiducial central model. The parame- 
ters and fiducial values we have chosen are described in the 
previous section. 

• A set of n measurements of the data x = {xi,X2, ... , Xn) 
(for e.g., the spherical hajmonic galaxy-galaxy power spec- 
trum C{£) over a range I = 1 ... n), and a model for how the 
data depend on cosmological parameters, i.e.: x — x (O ^ 

• An estimate of the uncertainty on the data A{x), which 
may depend on the given experiment (instrument noise, shot 
noise, etc..) as well on the data estimator (e.g., cosmic vari- 
ance) . 



The Fisher Information Matrix (FIM) is defined by: 



(Al) 



where C = — InL and L — L{x, 9) is the likelihood function 
or the probability distribution of the data x, which depends 
on some model parameter set 9. 

The forecasted uncertainty on the parameter 9i can be 
estimated directly from the FIM and obeys: 



(A2) 



By assuming the errors on the estimator x axe Gaussian 
and independent, then the Fisher matrix can be estimated 
by: 



51 (Aa;„)2 89, 



dx„ dx„ 



(A3) 



The above equation holds for a single redshift; the ex- 
pressions for the tomographic calculation are given below 
for each BAO method. 



Al Full Fourier space power spectrum P{k) 

For the case of the full Fourier space power spec- 
trum P(k), the Fisher matrix can be approximated a s 
(|Tegmark. Tavlor fc Heavenslll997l : ISeo fc Eisensteinll2003l ): 



It 



ainP(fc) 9lnP(k) -,rfj:^ dk 



2{27r)-l 



(A4) 



where the effective volume is given by: 
n{r)P{k,^^) 



Kff(fc,M) 



n{f)P{k,y) + 1 



(A6) 



and the matter power spectrum is calcu lated using the 
CMBFast code (jZaldarriaea fc Seliak|[l999h . 

Uncertainty is the redshift determination will leave tan- 
gential modes unaffected and s mear ra dial modes. This loss 
of information is modelled by (|Seo fc Eisenstein., 2003 ') : 



P(fc) = Pobs(fc)e" 



(A7) 



i.e. by convolving the radial position by a Gaussian uncer- 
tainty. 
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A2 Fourier space BAO 'Wiggles Only' 

In the case of the 'wiggles only' method, only the size of the 
radial and tangential BAO scale is used. This method omits 
all information from the shape and amplitude of the power 
spectrum, the redshift distortion and even the amplitude of 

the BAO wiggles. 

We follow the method o f ISeo fc Eisensteinl(|200'it . which 
we briefly summarise here. The starting point is the expres- 
sion to compute the Fisher matrix from the full Fourier space 
power spectrum f equation IA5|l . The non-linear power spec- 
trum used is extracted from the linear power spectrum with 
a damping term representing the loss of information coming 
from the non-linear evolutionary behaviour. This damping 
term is approximated by an ex ponential computed usin g the 
Lagrangian displacement fields lEisenstein et al.l (|2006h : 



Pni(fc, ^^) = -Pim(fc, m) exp( 



(A8) 



This non-linear damping can be decomposed into the 
radial and perpendicular directions by: 



Pnl(fc, y) = Plin(fc, m) exp 



(A9) 



A delta function baryonic peak in the correlation func- 
tion at the sound horizon scale, s, translates into a 'wiggles 
only' power spectrum of the form 



Pb oc 



sin(fcs) 
ks 



(AlO) 



This functional form is also obtained if the power spec- 
trum is divided by a smooth version of it or if we only use 
the b aryonic part of the pow er spectrum transfer function 
(e.g.. lEisenstein fc Hul l|l999l )'). However, the baryonic peak 
is widened due to Silk damping and non- linear effects. This 
Gaussian broadening translates into an exponential decay- 
ing factor in Fourier space: 



\{ks) 



ks 
sin(fcs) 



ks 



exp 



exp 



-(fcEsiik)^'*] exp 



-k^T.1 



(All) 
(A12) 



We can further refine the calculation to 2 dimensions, 
separating the radial and tangential location of the baryon 
acoustic pe ak in the correlati on f unction. Using the sam e 
notation as iBlake et al.) l|2006l) and |Parkins"on et al. I l|200'if ). 
the observables are y{z) = r{z)/s and y'{z) = r'{z)/s re- 
spectively, where r{z) is the comoving distance to redshift 
z. Measuring the fractional errors on these two observables 
is equivalent to measuring the fractional errors on H ■ s and 
Da/s. 

We can now use the derivatives of t he 'wiggles only' 

J iower spectrum by these quantities (see ISeo fc EisensteinI 
2003) for all the details) to compute the Fisher matrix 
(equation 26 of^S eo fc EisensteinI l|2007l )). This equation in- 
cludes the degradation due to redshift distortions in the fac- 
tor R{fi) = (1 -I- /3/i^)^ - though does not include any infor- 
mation from these distortions (see Table [T]). 

The effect of photometric redshift errors can be included 



in this formalism as an additional exponential term in the 
redshift distortion factor, R{fi) = (H~/9/i^)^ exp(— fc^/i^E^), 
where Ez is the uncertainty in the determination of photo- 
metric redshifts. We evaluate these integrals to obtain the 
Fisher matrix for the angular diameter distance, Da/s, and 
the rate of expansion, H ■ s. 

To compute the errors on the parameters, we follow 
IParkinson et all l|2007f ). with some changes: Since we know 
the correlation between the errors on the vector yi — 
{r{zi) /s,r' (zi)/ s} , we form a small 2x2 covariance matrix 
for each redshift bin, which we inv ert to obtain the corre- 
sponding Fisher matrix, F'''. While IParkinson et al.l (|2007l ) 
compute the transformation to the cosmological parameters 
analytically, we perform it numerically. 

The above is enough to provide constraints on cosmo- 
logical parameters. We also use another equivalent imple- 
mentation which adds another step and is described here- 
after: We form a Gaussian likelihood L oc exp(— x'^/2) with 



i a, 6—1 



(A13) 



where i runs over the redshift bins and a, b over r/s and r'/s, 
and A = Xdata(2i)~-'^theory(2i) is the difference between the 
data and the theory. We then compute numerically the ma- 
trix of second derivatives of at the peak of the likelihood, 
via finite difi'erencing in all parameters (where it is necessary 
to divide by 2 when using rather than — InL). 



A3 Spherical Harmonic space C{1) 

In spherical harmonic space the estimator is the galaxy- 
galaxy angular 2-point function C(i). The Fisher matrix is 
then calculated using: 



where 
ACs3(£) 



(2^+l)/,ky 



(A14) 



(A15) 



and A/'g is the galaxy shot noise. 

The galaxy 2-point correlation function in spherical har- 
monic space is given by: 

C,,{£) = {\ai^f) = W(zO y dk^^lWKk)^ (A16) 

where airn represent the spherical harmonic coefficients and 
the galaxy bias is taken to be constant across the depth of 
each redshift bin. The power spectrum can be expressed as: 

A'(fc) = (A17) 

and is calculated using the publicly available code CAMB 
(jLewis et al.ll2000l ). 

The real space window function Wlik) is given by: 



Wl{k) = / dre(r)j,(kr)D(r). 



(A18) 
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The normalized galaxy distribution is denoted by 0(r), the 
term jiir) refers to the spherical Bessel function of order I 
and D{r) is the growth function. 

In redshift space (see sect ion [?!3)) the 2-point correlation 
function has an extra term: 



F, 



CMB 



C{t) = 47rr 



dk 



A'(k) 



w^(k) + /3w^(k)|^ 



(A19) 



where B = is the distortion term which modulated the 
amplitude of the distortion and 

1 



W^{k) = ^ I dr^j;(kr) 



(A20) 



The reader is referred to [ Fisher et alj (|1994l ) for a derivation. 

iPadmanabhan et al.l l|2006D showed that this could be 
re- written: 



Wik) = J dre(r) {AeUM " ^iU-2M " D,jf+2(kr)) , (A21) 
where: 



Ae = 
Be = 



(2£ + 3)(2^- 1) 

(2£+ 1)(2^ - 1) 

(£ + l)(^ + 2) 
(2^ + l)(2f + 3) 



(A23) 



(A24) 



The total galaxy window including redshift distortions 
can then be rewritten: 

H/;+"(fc) = W''{k)+pW%k) , (A25) 
= W^ik) + 13 {AeW[{k) + BeW[_2{k) + DiW[+^{k)) (A26) 

So that the 2-point function with redshift distortions: 



C{£) = 47rb^ 



dk^^\W[+' 

ft 



(A27) 



APPENDIX B: 
PRIORS 



CALCULATION OF PLANCK 



The constraints from the baryon wiggles (alone or including 
the full galaxy correlation) can be combined with constraints 
from other probes. In this paper we exploit measurements of 
the CMB anisotropics as a complementary probe; combined 
with constraints from the galaxy survey, this will tighten 
constraints on some of the cosmological parameters. 

The primary constraint on cosmology from the CMB 
comes from the measurement of the angular size of the 
sound horizon at last scattering. We will use the forthcom- 
ing Planck mission as benchmark for a CMB prior. In order 
to forecast the ability of Planck to constrain cosmological 
parameters we will need to estimate the errors on the mea- 
surements of the temperature and polarization power spec- 
tra. We will conservatively not include any B-modes in our 
forecasts and assume in our fiducial model no tensor mode 
contribution to the power spectra. 

The Fisher matrix for CMB power spectrum is given by 
l|Zaldarriaga fc Seliak|[l997l : IZaldarriaga et al.lll997l '): 



EE 



(Bl) 



where 6i are the parameters to constrain, Cx,i is the har- 
monic power spectrum for the temperature-temperature 
(X = rr), temperature- E-polarization {X = TE) and the 
E-polarization-E-polarization {X = EE) power spectrum. 
The covariance COVxy of the errors for the various power 
spectra is given by the fourth moment of the distribution, 
which under Gaussian assumptions is entirely given in terms 
of the Cx I with 



h{CE,i + Wp^Biy 
fi Cte.i + 
{Ct,i + Wr^B-^) [Ce.i + W^p'^r") 

ftCTE,l 

fiCTE,i {Ct^ + Wt'b;-^) 



(B2) 
(B3) 
(B4) 



= fiCTE.l {Ce,i + Wp^Bf^) , 

and Wt,p = (crT,pSfwhi, 



(B5) 
(B6) 
(B7) 



COVt.t 
COVe.e 

,te 

COVt,e 
COVt ,TE 

,TE 

where ft = {2i+i)f^i^y ^'i'^ Wt,p = (crT,pSfwhm)~'^ is the 
weight per solid angle for temperature and polarization, with 
a Icr sensitivity per pixel of or.p with a beam of Sfwhm ex- 
tend. The beam window function is given in terms of the 
full width half maximum (fwhm) beam width by Be — 
exp (^—£(fi + l)9f„Y^^/16ln2^ and /sky is the sky fraction. 
Note that equation IBll usuallv includes a summation over 
the Planck frequency channels. However we conservatively 
assume that we will only use the 143 GHz channel as science 
channel, with the other frequencies used for foreground re- 
moval, which is not treated in this paper. This channel has 
a beam of Sfwhm = 7. 1' and sensitivities of or = 2.2iiK /K 
and (TP = 4:.2^iK/K ijThe Planck CoUaborat"ior] l2006h . To 
account for Galactic obstruction, we take /aky = 0.80. Note 
we use as a minimum ^-mode, £min = 30 in order to avoid 
problems with polarization foregrounds and subtleties for 
the modelling of the integrated Sachs-Wol fe effect, which 
depends o n the specific d ark energy model (|Weller fc Lewis! 
[2003; Caldwell fc Doranl i2005). 

We have now all the ingredients to calculate the Fisher 
matrix forecast for Planck CMB observations. However, we 
still have to specify the most suitable parameter set for do- 
ing this. It is well known that one of the primary parameters 
which can be constrained by the CMB anisotropies is t he an- 
gular size of the sound horizon (|Kosowskv et al.ll2002l ) . Also 
we should keep in mind that the Fisher matrix approach 
is a Gaussian approximation of the true underlying likeli- 
hood. In this context we would like to choose parameters 
in which the likelihood is as similar as possible to a Gaus- 
sian. Since the CMB anisotropies, with the exception of the 
integrated Sachs- Wolfe effect, are not able to constrain the 
equation of state of dark energy it is likely that if we would 
choose {wQ,Wa) as parameters for the CMB Fisher matrix 
this degeneracy is artificially broken by the Fisher matrix 
approach. This was also recognized by the DETF and we 
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hence f ollow their approach to calculate the CMB Fisher 
matrix l|Albrecht et al.ll2006l ). 

We choose as fiducial parameter set 9 = 
(wmjSs, In y4s,t<J6,ns,r), where 6s is t he angular size 
of th e sound horizon at last scattering (|Kosowskv et alj 
|2002| ). In As is the logarithm of the primordial amplitude 
of scalar perturbations and r is the optical depth due to 
reionization. Note that there might be an e ven more suitable 
param eter set, which is currently explored (|Mukheriee et all 
|2008| ). The optical depth in terms of the BAO analysis pre- 
sented here is a nuisance parameter and we marginalise over 
it analytically. We then calculate the Planck CMB Fisher 
matrix with the h elp of the publicly available CAMlfl code 
(|Lewis et al.ll2000l ). As a final step we transform the Planck 
Fisher matrix in the DETF parameter set to the one used 
in the analysis presented in this paper. The parameters 
here are given by ^ = (fi™, n^e, ft, crs, f^b, wo, Wa, ns) by 
using the transformation with the Jacobian 

J&a ~ — — (B8) 

and the Fisher matrix in our set of basis parameters 

(B9) 



This paper has been typeset from a TJjX/ I^Tj^X file prepared 
by the author. 



* http://camb.info 
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